Signal of Bose condensation in an optical lattice at finite temperature 
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We discuss the experimental signal for the Bose condensation of cold atoms in an optical lattice 
at finite temperature. Instead of using the visibility of the interference pattern via the time-of-flight 
imaging, we show that the momentum space density profile in the first Brillouin zone, in particular 
its bimodal distribution, provides an unambiguous signal for the Bose condensation. We confirm 
this point with detailed calculations of the change in the atomic momentum distribution across 
the condensation phase transition, taking into account both the global trapping potential and the 
atomic interaction effects. 
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A recent work raised an interesting question on how 
to unambiguously confirm the Bose condensation asso- 
ciated with a superfluid phase at finite temperature for 
cold atoms in an optical lattice [1, 2]. The signal of Bose 
condensation in an optical lattice is typically connected 
with the interference peaks in the time-of-flight images. 
It is shown in [1] that for a thermal lattice gas above the 
condensation temperature, one could also see the interfer- 
ence peaks with a pretty good visibility [3] . This thermal 
visibility goes down considerably if one takes into account 
a number of practical effects, as shown in the recent work 
0]. On the other hand, however, the thermal visibility 
above the transition temperature is not negligible, and 
its value depends on the details of a number of system 
parameters. There is no explicit critical value yet for the 
interference visibility to separate the condensation region 
from the thermal region. 

In this work, we suggest to use the momentum space 
density profile in the first Brillouin zone (measured from 
the time-of-flight imaging) as an alternative method to 
characterize the Bose condensate in an optical lattice. In 
particular, similar to the free-space case, the bimodal dis- 
tribution of the atomic momentum distribution should 
provide an unambiguous signal for the Bose condensa- 
tion. To show the practicality of this method, we calcu- 
late the finite-temperature atomic momentum distribu- 
tion in a three-dimensional (3D) optical lattice, confirm- 
ing a distinctive change in the distribution across the 
condensation phase transition. There are several com- 
plexities in this calculation: first of all, one needs to take 
into account the inhomogeneity of the lattice due to the 
global harmonic trap, which tends to broaden the mo- 
mentum distribution of the gas, in particular for the con- 
densate part. Secondly, the atomic interaction leads to 
two competing effects: it broadens the real-space distri- 
bution and thus sharpens the momentum space density 
profile in a trap; while at the same time, it also leads 
to the broadening of the momentum distribution during 



the time-of-flight process. In our calculation, we take all 
these effects into account. 

To illustrate the basic character of the experimental 
signal, first we calculate the finite temperature momen- 
tum distribution for free bosons in an optical lattice with 
a global harmonic trap, in which case exact solutions are 
possible after introducing some tricks. With the Fesh- 
bach resonance technique, one can also directly test the 
predictions in this case by turning off the atomic inter- 
action. Then we use the Hatree-Fock-Bogliubov-Popov 
(HFBP) approximation [1, 0, 0] to deal with the interac- 
tion between the atoms. The HFBP method is expected 
to be a good approximation for weakly interacting bosons 
except for a small region across the transition from the 
superfluid phase to the normal phase [5] . For cold atoms 
in a lattice with a global harmonic trap, such a region 
corresponds to only a slim layer, which has negligible in- 
fluence on the integrated momentum distribution. So we 
expect the HFBP approximation should be a reasonably 
good method for the calculation of the atomic momen- 
tum distribution before the Mott transition shows up in 
the trap. We note that the HFBP approximation has 
been used before for cold atoms, with the focus on the 
homogeneous case [6] or the one-dimensional system 

For free bosonic atoms in an inhomogeneous opti- 
cal lattice with a global harmonic trap, the Hamil- 
tonian has the form H = J dT'^\T)[—{h'^ /2m)d^ + 
Vop + F(r)]^(r), where m is the atomic mass, Vop = 
^ 2 3 ^^^^('^^*/^) optical lattice potential 

with d as the lattice spacing, and V^(r) = ^mcj^r^ is the 
global harmonic trapping potential. To diagonalize the 
Hamiltonian, we expand the atomic field operator ^(r) 
as ^(r) = ^pj^ w{r — R)aR, where the operator an {an) 
annihilates (creates) a particle on site w{r — R) is the 
Wannier function on the site R associated with the lat- 
tice potential Vop^ and the summation runs over all the 
lattice sites. A Fourier transform of this expansion gives 



^(k) =^(k)ak, 



(1) 
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where ^(k), iL'(k), ak denote respectively the Fourier 
components of ^(r), w{y)^ aR in the momentum space. 
Under typical experimental situations, the global har- 
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monic trap V{r) varies slowly on individual lattice sites. 
In that case, under the expansion above, the Hamiltonian 
is recast into the form 



k R 



(2) 



where the summation over k runs over the first Bril- 
louin zone. We have assumed here that the system 
temperature is well below the band gap, so the atoms 
only occupy the lowest band, where is well approx- 
imated by = — 2t ^^^-L 2 3 cos(/ci(i). The parameter 
t denotes the tunneling rate over the neighboring sites, 
with t ^ {3.5/y/7T)VQ^^e-'^^ for an optical lattice 0, 
where both t and Vq are in the unit of the recoil energy 
Er = h^7r^/2md^. 

The Hamiltonian (2) in principle can be directly diag- 
onalized numerically. However, the calculation of the fi- 
nite temperature momentum distribution requires deter- 
mination of all the eigenstates of the Hamiltonian, which 
is very time consuming for a three-dimensional lattice 
with many sites. Here, we adopt an approach which al- 
lows easy calculation of any finite temperature momen- 
tum distributions [1, [13, El- Note that the indices R 
and k in Eq. (2) are reminiscent of the coordinate and 
the momentum variables in quantum mechanics. So we 
can write down a first quantization Hamiltonian corre- 
sponding to Eq. (2) in the momentum space, where R is 
replaced by the momentum derivative d\^. The resulting 
Hamiltonian takes the form 



(3) 



This Hamiltonian describes free particles with effective 
mass m* = fi^/{muo'^) trapped in an effective poten- 
tial E\^ with periodic boundary condition (the period is 
given by the reciprocal lattice vector K). The quasi- 

momentum distribution (^(i\^ci\^ is then given by the 

square of the eigenstate wavefunctions (k)|^ of i^e//, 
averaged over all the eigen-levels n with a Bose distri- 
bution factor ^(en) = 1/ exp [(en — /i) /T — 1] at finite 
temperature T, where en is the corresponding eigenen- 
ergy and ji is the chemical potential fixed by the total 
number of atoms. From Eq. (1), we then obtain the 
atomic momentum distribution 



^(k) = <^t(k)^(k)) = mk)|^(atak; 

;(k)|^^^(en)|(/>„(k)|^ (4) 



The signal from the time-of-flight imaging corre- 
sponds to the column integrated momentum distribution 
'n±{kx^ky) = / n{k)dkz. 

Through Eq. (4), we have calculated the momentum 
distribution of a free Bose gas in an optical lattice with 
a global harmonic trap at different temperatures around 



the condensation transition. The parameters are chosen 
to be close to those of a typical experiment for ^''Rb, 
with the lattice barrier Vq = I^Er and the total particle 
number N :^ 10^. The results for the integrated mo- 
mentum distribution n±{kx, ky) are shown in Fig. 1 un- 
der different trap frequencies. In a weak harmonic trap, 
there are distinctive interference peaks even for a ther- 
mal gas above the transition temperature (Fig. 1(a)), 
which agrees with the result of Ref. [l| for a homoge- 
neous optical lattice that corresponds to the zero-trap 
limit. These peaks are caused by the short range ther- 
mal correlations between different sites. The correlation 
function of the thermal gas decays exponentially with 
distance x by the form e~^/^, where the characteristic 
length L is estimated to be ~ 1. Id for Fig. 1(a). As the 
trap frequency increases, the thermal interference peaks 
become less visible and eventually disappear (Fig. 1(c)), 
which agrees with the recent calculation of the visibil- 
ity in Ref. [J]. In the case of a weak global trap, it 
gets difficult to use the visibility to distinguish the con- 
densation phase transition as this parameter has a pretty 
high value for both the condensed and the non-condensed 
phases near the critical temperature [H, I3]- However, if 
one directly compares the interference patterns in Figs. 
1 (a,c) and (b,d) across the condensation transition, the 
difference in the interference peaks is still pretty clear: in 
particular, for the case with a condensate component, the 
interference peaks become much sharper. This reminds 
us to look at the momentum space density profile in the 
first Brillouin zone, which, compared with the single pa- 
rameter of visibility, should be able to give more detailed 
information about the system and its phase transition. 

With this in mind, we plot the column integrated 
momentum distribution along one of the axes passing 
through the center of the first Brillouin zone. The results 
are plotted in Fig. 2 (a) (c) (e) for different condensate 
fractions. It is clear that while the density profile is a 
thermal distribution for a normal e: as, a bimodal struc- 
ture appears when a non-zero condensate fraction exists 
in the lattice. Based on this observation, we propose to 
use the interference peaks as well as the bimodal struc- 
ture in the first Brillouin zone as an unambiguous signa- 
ture for detection of the Bose condensation in an optical 
lattice. 

The calculation above is useful for an illustration of 
the general qualitative features of the system, but we ne- 
glect the interaction between the atoms in the lattice. 
In a more realistic model, such interactions need to be 
taken into account, as they can greatly modify the time- 
of-flight images. Firstly, the repulsive interaction will 
tend to broaden the spatial distribution of the atoms and 
hence narrow the momentum distribution in the trap. 
Secondly, in the process of the time-of-flight expansion, 
the remnant atomic interaction transforms the interac- 
tion energy into the kinetic energy, which shows up in 
the final imaging process and broadens the momentum 
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FIG. 1: The intensity plot for the column integrated mo- 
mentum distribution near the first Brillouin zone with differ- 
ent trapping frequency a;/27r and temperature T. Figs, (a) 
and (c) are for the normal phase, and (b) and (d) are for 
the superfluid phase with the corresponding parameters: (b) 
T/Tc - 0.7; (d) T/Tc - 0.94. The total particle number is 
fixed to be Nt = 10^ for all the cases, and the lattice barrier 
Vb = IOEr. The unit for the axes is 27v/d. 



distribution. In the presence of interaction, the system 
can not be solved exactly. In the following, we calculate 
the momentum distribution of an interacting Bose gas in 
an optical lattice with a global harmonic potential trap 
using the Hatree-Fock-Bogliubov-Popov (HFBP) approx- 
imation [5]. 

To calculate the momentum space density profile for an 
atomic gas in an inhomogeneous optical lattice, we com- 
bine the HFBP method with the local density approx- 
imation (LDA). Under the LDA, a local region of the 
harmonic trap is treated as a homogeneous lattice sys- 
tem. We first look at the Hamiltonian with inter-atomic 
interaction for a homogeneous lattice gas: 



H 



-t ttj ttj + — ajttj ttitti — /i ajtti (5) 

(iJ) i i 



where t is the tunneling rate defined before, /i is the lo- 
cal chemical potential, U = U^g J \ w{r)\^dr is the on-site 
interaction rate with an approximate form of U/ Er ^ 
3.05{Vo/ER)^-^^{as/d) (Uhg is related to the s-wave scat- 
tering length by U^g = 47rfi^as/m, and = 5.45nm for 
^^Rb [H, I3l). We then transform all the field operators to 



the momentum space, and write the momentum compo- 
nent ttk as ttk = v5\^o + Jttk with the standard HFBP ap- 
proach, where v = (a^) = (ao) represents the condensate 
fraction, and 5a\^ and Jaj^ represent the excitations above 
the mean field. We keep the Hamiltonian to the quadratic 
order of the operators 6a\^ and Jaj^, and then diagonal- 
ize it to get the thermal potential Q = —T\ntr (e~^/-^). 
The stationary condition dVt/dv = gives the saddle 
point equation for the chemical potential: 



fi = Eo — Uno ^ 2Un, 



(6) 



where = —^t for a three-dimensional system, no is 
the per-site density of the condensate fraction, and n is 
the total number of particles per site. With the thermal 
potential, it is easy to derive the number equation and the 
atomic momentum distribution. For the non-condensate 
part with k 7^ 0, the per-site quasi- momentum density 

distribution (aj^ak) is given by 



aj^ttk 



- £^0 + Unp ^^^^ f hcj^ 



where hij\^ = ^/ {E\^ — Eq)'^ + 2Unp{E\^ - Eq) is the dis- 
persion relation for the Bogliubov excitations. With 
the LDA, the local chemical potential /i(r) at a dis- 
placement r from the trap center is determined from 
/i(r) = jiQ — F(r), where /io is the chemical potential 
at the trap center. We fix the total number density per- 
site at the trap center, which, together with Eqs. (6-7) 

and the number equation n = no + Xlk/o ^^k^k^, al- 
low us to calculate the chemical potential /io- We may 
then determine /i(r) and ^<^k«k^ at any trap location r. 
The overall non-condensate part of the atomic momen- 
tum distribution is given by the integration of ^<^k<^k^ 
over the whole harmonic trap, multiplied by the Wan- 
nier function |K;(k)|^ as shown in Eq. (4). 

For the condensate part, the LDA will lead to an ar- 
tificial ^-function at zero momentum. To avoid this ar- 
tifact, one needs to consider explicitly the broadening 
of the condensate momentum distribution by the har- 
monic trap. From the above LDA formalism, we get the 
condensate fraction no(r) at any trap location r. The 
condensate wavefunction in the trap can be well approx- 
imated by (j)o{v) = Y^no(r) (which actually corresponds 
to the solution of the Gross-Pitaevskii equation under 
the Thomas- Fermi approximation [lH). The condensate 
part of the atomic momentum distribution is thus given 
by the Fourier-transform of the wave function 0o (r) . The 
resulting momentum distribution is then added to the 
momentum profile of the normal gas which is typically 
at the edge of the trap. We plot the results of our calcu- 
lation in Fig. 2 (b)(d)(f). In general, the peaks are now 
sharper compared to the case of an ideal Bose gas, which 
is just as expected. A bimodal distribution still shows 
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FIG. 2: (Color online) Column integrated momentum density 
profile (shown along the x-axis) for (a)(c)(e) non-interacting 
lattice bosons and (b)(d)(f) interacting lattice bosons. The 
solid curves are the total momentum density profile; the thick 
dashed curves are the momentum density profile of the nor- 
mal component; the dotted curves are the momentum den- 
sity profile taking into account of the interaction broaden- 
ing during the time-of-flight evolution, (a) No/Nt ^ 0.46, 
T/Tc ^ 0.75 (A^o is the total number of atoms in the conden- 
sate and Nt is the total number of atoms in the trap); (b) 
No/Nt - 0.52, T/Tc - 0.34; (c) No/Nt - 0.10, T/Tc - 0.95; 
(d) No/Nt - 0.10, T/Tc ^ 0.65, inset: enlarged density profile 
below the dotted line; (e) T - l.lTc; (f) T - l.lTc. For the 
non-interacting calculations, the total particle number is fixed 
to be Nt — 10^, and the trapping frequency is uj/2ii — 20Hz; 
for the calculations with interactions, the trapping potential 
is Lj/27T = 20Hz, the number density per-site at the center 
of the trap is fixed to be n(0) = 1, with the total number of 
particles on the order of Nt = 10^ ~ 10^. The lattice barrier 
in all the cases is Vb = IOEr. The unit for the x-axes is 27i/d. 



up when the condensate fraction becomes non-zero (Fig. 
2(d)). In real experiments, due to the finite resolution in 
imaging, the height of the central condnesate peak will 
be significantly reduced, and the bimodal structure will 
become more pronounced and should be readily observ- 
able. 

After turnoff of the trap, the atomic momentum distri- 



bution will be broadened during the time-of- flight process 
due to the collision interaction. The interaction broad- 
ening of the momentum distribution happens dominantly 
through the condensate fraction [3, Q , for which one has 
a higher number density and a smaller expansion speed. 
To incorporate this effect, we use the Gross-Pitaevskii 
equation [12] to solve the evolution of the condensate 
wavefunction 0o(r,t), starting at ^o(r) when the expan- 
sion time t = 0. The Fourier component of ^o(r,t) after a 
long enough expansion gives the final momentum distri- 
bution for the condensate part. The results of the calcu- 
lation are shown in Fig. 2 (b) (d) for ^^Rb atoms. One can 
see that the interaction broadening leads to some quanti- 
tative corrections of the profile by lowering its peak value, 
but it does not change much the overall picture. In par- 
ticular, the bimodal structure remains similar when there 
is a condensate fraction. 

In summary, we have shown through explicit calcula- 
tion that the interference peaks, combined with the bi- 
modal profile of the central peak in the first Brillouin 
zone, provide an unambiguous signal for the condensa- 
tion phase transition for cold atoms in an optical lattice 
at finite temperature. We develop some techniques to 
calculate the atomic density profile, in particular for an 
inhomogeneous system with a global harmonic trap. In 
general, the density profile of the interference peak gives 
more detailed information of the system, and a compari- 
son of the density profiles from the theoretical calculation 
and from the experimental observation will contribute to 
the understanding of this strongly correlated system at 
finite temperature. 
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